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Gravitomagnetic Field of Rotating Rings 

Matteo Luca Ruggiero^ 


Abstract In the framework of the so-called gravi- 
toelectromagnetic formalism, according to which the 
equations of the gravitational field can be written in 
analogy with classical electromagnetism, we study the 
gravitomagnetic field of a rotating ring, orbiting around 
a central body. We calculate the gravitomagnetic com¬ 
ponent of the held, both in the intermediate zone be¬ 
tween the ring and the central body, and far away from 
the ring and central body. We evaluate the impact of 
the gravitomagnetic held on the motion of test particles 
and, as an application, we study the possibility of using 
these results, together with the Solar System ephemeris, 
to infer information on the spin of ring-like structures. 


1 Introduction 

In General Relativity (GR) mass currents give rise to 
gravitomagnetic (GM) helds, in analogy with classical 
electromagnetism: actually, the held equations of GR, 
in linear post-newtonian approximation, can be written 
in form of Maxwell eq uations for the gravitoe l ectro- 
magnetic (GEM) held s ( Ruggiero and Tartaelia 2002 


Mashhoon et al.H2001al) . ( Mashhoon 2003 ). where the 
gravito-electric (GE) held is just the Newtonian held. 
Even though these effects are normally very small and 
hard to detect, there have been many efforts to mea¬ 
sure t hem. For instance, the fa mous Lense-Thirring 
effect ( Lense and Thirring IQlSl) . that is the preces¬ 
sions of the node and the periapsis of a satellite ornit- 
ing a central spinning mass, has been analyzed in dif¬ 
ferent conte xts: there are the LAGEOS tests around 
the Earth ( Ciufolini and Pavlid 2004 : Ciufolini et al. 
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2010a ). the MGS tests around Mars ( loriollioO^ , 2010 ) 
and ot her t ests around the S un and the planet s (llorid 
2012a ): see Ciufolini ( 2007 ); lorio et al. ( 2011 , l2013l) : 


Renzettil (|2013bf) for a discussion and a review of the 
recent results. In Fe bruary 2012 the LARES mission 
( Ciufolini et al. 112012 ) has been launched to measure the 
Lense-Thirring effect of the Earth, and is now gather¬ 
ing data; a com prehensive d i scussi o n on this mission 
can be found in loric ( 200^ 20091) : Renzetti ( 2013a . 
2OI2I) : Ciufolini et al. ( 2010bl. 2015 ). In th e recent past, 
the Gravity Probe B (jEveritt et al.ll201ll) mission was 
launche d to measure the prece ssion of orbiting gyro¬ 
scopes ( Pugh 1959 : Schifl 1960l) . The GM clock effect, 
that is the difference in the proper periods of standard 
clocks in prograde and retrograde circular orbits around 
a rot ating mass, has been inv e stigated but not detected 
yet 


( Mashhoon et al 

2OO2I: 


lorio et al 


2001 


Lichtenegger et al. 


Mashhoon et al 

200f)l) . 


2001b 


non¬ 


standard for m of g r avitomag netism has been recently 
analyzed by Acedol ( 2014al fh). in a purely phenomeno¬ 


logical context. Eventually, the possibility of testing 


sidered by many authors in the oast (Braginsky et al. 

1977, 

1984; Cerdonio et al. 1988i Liubicic and Logan 

1992; 

jCamacho 2001 

lorio 

2003 

: Pascual-Sanchez 

2003; 

IStedman et al. 

20031: 

lorio 

2006bl): a recent 


propo sal pertains to t he use of an array of ring 


lasersI Bosi et al. 2011 : Ruggierol 2015 ). and 


IS now 


underway (|Di Virgilio et al.ll2014ri. 

In a recent paper ( RuggierJ 2015 ), we have investi¬ 


gated the gravitational field of massive rings: exploit¬ 
ing the GEM analogy, we have studied both the GM 
and the GE components of the field, produced by a 
thin rotating ring, orbiting the central body along a 
Keplerian orbit. The ring field can be dealt with as a 
perturbation of the background field determined by the 
central body. We have used a power series expansion 
to calculate the field in the intermediate zone between 
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the central body and the ring. Massive rings are ubiq- 

equations take the form 

Giufolini and Wheeler 

1995; 

uitous and important in astrophysics, as suggested in 

Ohanian and Ruflini 

2013 



TorinI (I2m2hh : in iR.amos-Caro et aD (j201l[) the effects 


of geometrical deformations on ring-like structures are 
studied, together with the implications for stability and 
regularity of the motion of test particles (also for Sat¬ 
urn’s and Jupiter’s rings). Henc e, motivated by the 
relevance of ring-like structures, in Ruggierol ( 20151) we 


have focused on the GM compon ent of the hel d (the 
GE one is exhaustively studied in loriol ( 2012bl) ). and 


studied its impact on some gravitational effects, such as 
gyroscopes precession, Keplerian motion and time delay 
in some simplified geometric configurations. The under¬ 
lying idea is to consider the possibility of using these 
tests to estimate the mass and the angular momentum 
of matter rings. Here, we want to pursue the study of 
the GM field of rotating rings: to be specific, we want 
to calculate the GM field in the whole space, both in 
the intermediate zone between the ring and the central 
body and far away from the ring and central body. As 
for the effects of the ring field, we will focus on the per¬ 
turbations of the Keplerian orbital el ements of a test 
particle: while in the previous paper Ruggier3 ( 2015 ) 


we have considered just the case of coplanarity between 
the ring and the test particle orbit, here we will con¬ 
sider an arbitrary configuration. Then, we will compare 


vations of Solar System ephemeris (Fienga et al. 

2011 

Pitieva and Pitiev 

2013; 

Pitiev and Pitieva 2013 

). 


The paper is organized as follows: we review the 
foundations of the GEM formalism in Section[21 while in 
Section |3] we obtain the GM field of the ring; in Section 
[4] we focus on the perturbations of the orbital elements 
determined by the GM field, and use the recent data of 
Solar System ephemeris to estimate the spin of ring-like 
structures. Gonclusions are eventually in Section [5] 


2 The GEM formalism 

If we work in the weak-field and slow-motion ap¬ 
proximation, we may write the space-time metric 
in the fornj^ in terms of the 

Minkowski tensor and the gravitational poten¬ 
tials which are supposed to be a small pertur¬ 
bation of the flat space-time metric: -C \g^v\- 

Hence, in linear approximation, on setting = 

hfiu — with h = tr(/i^y), and imposing the 

transverse gauge condition = 0, the Einstein 


^Greek indices run to 0 to 3, while Latin indices run from 1 to 3; 
bold face letters like x refer to space vectors. 


^ - IGttG^ 

(—1 tl^i^ = ^ 1 . 


( 1 ) 


It is a well known fact that, due to the a n alogy with 
electromagnetism (MashhoonI 200lL 2003t Mashhoon 

^ ' iKnmafl^^ 


200, 5 I 1199,4 iMashhoon et al 


et 


a1.ll200Sh 


the solution of the field equations © can be written in 
the fornll 


ds^ = —^1 — 2-^^ dt^-(A • dr)dt-I- 


(^1 + 2 ^^ Sijdx'-dx^ 


( 2 ) 


in terms of the gravitoelectric ip (/iqo = 4-^ ) and grav¬ 
ed 

itomagnetic Ai [h^i = —2^) potentials, which are re¬ 
lated to the sources of the gravitational field by 


p{ct,v) = / 
Jv 


p{ct — |r — R|, R) 


dH, 


c Jv 


Ir-RI 


dH. 


(3) 


(4) 


In the above equations p is the mass density and 
is the mass current of the sources. So, we see that, 
besides the usual Newtonian contribution p, related to 
the mass of sources, there is a contribution related to 
the mass current of the sources. The gravitoelectric E 
and gravitomagnetic B fields are then defined as 


1 5A _ 


B = VA A 


(5) 


For stationary sources, the equation of motion (i.e. the 
spatial components of the geodesics) of a test mass 
mtest moving with speed v in GEM fields E, B turns 
out to be (see e.g. Bini et all ( 20081) ') 


dv 


f^test j, — ^testE 2?7itest 

dt C 


- X B, 


( 6 ) 


to lowest order in v/c. In the convention used, a 
test particle of inertial mass mtest has gravito-electric 
charge qe = —mtest and gravito-magnetic charge qb = 
—2mtest; the GEM Lorentz acceleration acting on a test 
particle is 


A = -E - 2- X B 


(7) 


^Here and henceforth we u se the convention introduced by 
Mashhoon iMashhoonl Il993l ) to exploit the standard results 
of electrodynamics to describe gravity in post-Newtonian lin- 
ear approximation. Othe r conventions are used e lsewhere 
jCiufolini and Wheeleilll995l : lOhanian and RufHnil 1^131 ). 
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3 The gravitomagnetic field of rotating rings 



Fig. 1 The point P has spherical coordinates (r, 8, (j>), and 
the origin of the coordinate system coincides with the cen¬ 
tral body; they are related to the Cartesian ones by r = 
x'^ + + z'^ , (f) = arctan(y/a;), 9 — arctan(^a;^ -I- y^jz). 

The ring is in the xy plane, where 9 = 7r/2. 


In this Section, we calculate the gravitomagnetic 
field produced by a rotating ring; we suppose that the 
ring is thin and made of continuously distributed mat¬ 
ter with constant density, orbiting a central body. Fur¬ 
thermore, for the sake of simplicity, we assume that 
the ring is circular: actually, the case of an ellipt ically 
shaped ring has been considered in iRuggierol (|2015l) , but 
the resulting expressions are in general unmanageable, 
even to lowest order in the eccentricity. 

The central body is supposed to produce its gravita¬ 
tional field, which is determined by its mass M, and 
angular momentum S. In the inertial frame where 
the central body is at rest, we set a Cartesian coor¬ 
dinate system {x,y,z}, with the corresponding unit 
vectors Ua;,Uy,u^; if the body is located at the ori¬ 
gin and its angular momentum is directed along the 
z axis, S = Su^, the space-time metric to lower order 

has the form (ED, with (f =- 

r 

= \r\ = 


A G(SAr) 

A =-=—, where 


r = |r| = va; 

turns out to be 


+ y'^ + z^. In particular, the GM field 


B 


G 

c 


3(S • r)r 


( 8 ) 


and has a typical dipolc-like behaviour: in other words, 
it is analogous to the magnetic field produced by a 
dipole. 

In order to evaluate the gravitomagnetic field of the 
rotating ring, we proceed as follows: an infinitesimal 


mass element dm of the ring is orbiting around the 
central body; we know the total mass m of the ring, 
and its angular momentum s, which we assume to be 
constant: in other words, we consider a stationary 
ring. In our perturbative approach, we do suppose that 
m <g; M, s <C 5*. Due to the presence of this ring, 
the gravitomagnetic potential (|1]) is perturbed, so that 
A —>■ A -I- A''™s^ where |A’''"®| <C |A|. In particular, 
we are interested in calculating this perturbation, by 
means of a power series expansion, (i) in the interme¬ 
diate region between the central body and the ring, (ii) 
in the outer region of the system, i.e. away from the 
ring and the central body. 

To this end, we consider the following geometric con¬ 
figuration: we suppose that the ring is in the xy plane, 
that is the symmetry plane of the central body. In order 
to deal with the symmetries of the problem in a sim¬ 
pler way, we will also use spherical coordinates {r, 0, ^}, 
togheter with the corresponding unit vectors u^, ug, u^. 

Let P denote the point where we want to evaluate the 
GM field (see Figured]): its spherical coordinates are 
{r, 6, (f)} and its position vector is r; the position vector 
of a mass element dm of the ring is R, where R = |R| 
is the radius of the ring, and its spherical coordinates 

771 

are {R, 7r/2, dJl.The ring uniform density is A = -——; 

2t:R 

furthermore, v is the (constant) modulus of the mass 
elements speed. We must substitute fdV —>■ AtiML in 
Q: P are the components of the velocity, which may 
be written as v = — nsin^Uj; + 'ccos'hUj,, and dL is 
the infinitesimal arc length of the ring. Accordingly, we 
get 


A''™® = 


2G 


A (—n sin $ Uj; -I- n cos $ Uj^) dL 


(9) 


We may write dL = Rdd>; hence, on introducing the 
angular momentum per unit mass a = v R, the above 
integral can be written as 


A’''"® = 


2GAct ^’^ (—sin® Uj,-I-cos $ Uj^) d$ 

c y$=o Ir - Rl 


( 10 ) 


This expression can be expanded in power series: for 

T 

(i) r < R, we expand in powers of e = while for(ii) 

R 


r > R, we expand in powers of e = —. Consequently, 

r 

we may write |r — R| in the form 


{i)r < R: 
(ii) r > R : 


|r — R| = R)J 1 — 2^ sin0cos((/) — $) -|- 


1 — 2— sin0cos(()) — $) -f 
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or 

|r — R| = d\/l — 2e sin 0 cos{4> — <&) + (11) 

, {i)r < R : d = R, e = r/R 
where: 

[ii)r > R : d = r, e = R/r 
Hence, we may write 

2GXa 1''^='^^ (—sin$ u^; + cos<I> Uy) d<i> 
cd Y^l — 2e sin 0 cos{(j) — $) + 

( 12 ) 

Because of the cylindrical symmetry, we may choose 
the observation point at ^ = 0: as a consequence the 
X component of is null, and the y component is 

equal to 


According to Eq. ®, the corresponding gravito- 
magnetic field can be obtained from = V A A“'™s. 
For r <C i?, the gravitomagnetic field has the following 
components 


^ yCrS ^ ^ r 

^ cos 9 -cos 0 sm0— + 

cR-^ cR'^ R 


-f 

Gs 

—^ cos 0 (30 sin^ 0 — 3 



)(16) 

= 

2Gs . ^ 9Gs . 2 

r 

e-^ 

a 



-f 

Gs 

—^ sin 0 (6 — 30 sin^ 0 

)^ + o( 


1(17) 


While, for r ^ R, we obtain 


= 


2GXa 


l•^=2^v 


COS $ d$ 


(13) 


cd 7$=o y^l — 2esin0cos(<i)) + 

As shown in .Tacksonl (199^, the above integral (|T^ 


can be evaluated in terms of elliptic integrals: 


f.<J>=27r 


COS $ d$ 


&=o ^1 — 2e sin 0 cos(<I)) + ( 


VTT^ 


2e sin 0 
4e sin 0 


{2-p^)K{p)-2E{p) 


(14) 


where --—;—r and K(p), E(p) are the 

complete elliptic integrals of first and second kind. For 
e <C 1 (that is for i? <C r or r <C R) the result of 

1 • 1 • jiTT 1 7resin0 

the integral m Fq. (1141) becomes - tttt 

(l + e2 + 2esin0)^/" 
and, consequently, we may write the gravitomagnetic 
potential in the form 


AT^ = 


2GXa 


TTC Sin t 


cd (1 + e2 + 2esin0)^^^ 

If we perform a power-series expansion we obtain: 


(15) 


(*)r«R:A™8 = 


Gs sin 0 r ^ Gs sin^ 0 


-f 


cR^ R 
15 Gs sin^ 0 

y 


ci?2 


ci?2 R2 

3 Gs sin 0 
2 


ci?2 


r 


{ii)T^R: A^l^^ = 




Gs sin 0 „ Gs sin^ 0 R 

2 ^ 2 

cr or r 

,3 


+ O 


15 Gs sin'^ 0 
2 cr'^ 


3 Gs sin 0 

2 cr2 


R^ 


= 


2Gs 9Gs R 

— 5 - COS 0 5 - cos 0 sm 0— 


cr"’ 


cr"” 


cos 0 (30 sin2 0 — 3) ^ + O (^ 


(18) 


= 


Gs 6Gs •} „R 

■ sm0- 5 - sm 0- 

Gs 


i?2 






We notice that, to lowest approximation order, the 
gravitomagnetic field has the following expressions 


(*)r<i?:B™s = 


{ii) r i? : B'' 


2Gs 

ci?3 

2G 

ci?3‘ 

CJ.3 

G 
c 


(cos0Ur — sin0ue) = 


( 20 ) 


(2 cos 0Ur -f sin Bug ) = 


3(s • r)r 


( 21 ) 


The expression (1^01) of the field i nside the ring i s in 
agreement with the one obtained in Ruggierol ( 2015ll : on 
the other hand, we see that the expression (EH of the 
field outside the ring is, as expected, the usual dipole 
field. In the following Section we are going to use these 
expressions to calculate the perturbing acceleration on 
the motion of orbiting bodies and, then, the correspond¬ 
ing variations of the orbital elements. 


4 Gravitomagnetic perturbations 

In this Section we evaluate the impact of the GM of the 
rings on the orbital elements of test particles: in par¬ 
ticular, we consider below the effects on Solar System 
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bodies. To this end, we use the expressions of the GM 
held (l20l) and (|2T|) to calculate the perturbing acceler¬ 
ation 

W = -2- X (22) 

c 

then, we can evaluate its effects on planetary motions 
using the Gauss equations for the variations of the el¬ 
ements, which enable us to study the perturbations of 
the Keplerian orbital elements due to a generic perturb¬ 
ing acceleration. 


said before, the ring is assumed to be stationary: this 
amounts to saying that the motion of the ring matter 
is constant during the particle’s timescale. 

We obtain the following results. 

For the test particles orbiting outside the ring, on 
using the expression (I 2 T]) of the GM held, we have non 
null secular variations only for the argument of perias- 
tron and the node: 


< d} >= — 


6Gs 1 
c2 a3(l-e2)^/^ 


(23) 



Fig. 2 Unperturbed orbit of the test particle 

To begin with, let us describe the conhguration of 
the unperturbed test particle orbit. We refer to Figure 
O besides the already mentioned Gartesian coordinate 
system {x,y,z}, we introduce another Gartesian coor¬ 
dinate system {X,Y,Z}, with the same origin. The 
ring lies in the xy plane, while the the orbital plane of 
the test particle XY plane. In particular, we denote 
with the angle between the x axis and the line of the 
nodes, while the angle between the z and Z axes is i. 
The periastron is along the X axis, and we denote by oj 
the argument of the periastron, i.e. the angle between 
the line of nodes and the X axis. 


orbita l elements (see e.g. 


Bertotti et al. ( 

2003 

), 

Roy 


verse (in-plane components) and normal (out-of-plane 
component) projections of the perturbing acceleration 
(j22p on the orthonormal frame comoving with the par¬ 
ticle, and then we use the Gauss equations for the vari¬ 
ations of the semi-major axis a, the eccentricity e, the 
inclination i, the longitude of the ascending node fl, 
the argume nt of peric entre oj and the mean anomaly 
M (see e.g. Rov ( 2005ll l. We want to stress that, as we 


<n>= 


2Gs 1 


(24) 


a3(l-e2)^/^ 

Eventually, for the longitude of the pericenter we have 
2Gs (3cosi — 1) 1 


< vu >= —- 


a3 (1 - e2)3/2 


(25) 


On the other hand, for test particles orbiting inside the 
ring, on using the expression (j20ll . we have the following 
non null secular variations: 


< O >= 


2Gs 

c2i?3 


< M >= n ■ 


2Gs cosi 


(4V^ 


— e2 -I-1 


(26) 


(27) 


Moreover, for the longitude of the pericenter we have 
2Gs 


< w >= 


2R3 


(28) 


We remember that fot an unperturbed Keplerian ellipse 
in the gravitational field of a body with mass M, it is 
n= v'GM/a3. 

The above results can be used to make a compar¬ 


ison with the recent ob s ervations (iFienga et al.l 2011 , 
Pitieva and Pitiev 2013 : Pitiev and PitievJ 201311 : f^ 


instance, on using the available supplementary ad¬ 
vances Azu, we may give estimates on the spin of ring¬ 
like structures in the Solar System. Let us start from 
planets orbiting outside the ring; in particular, we con¬ 
sider a hypothetical ring of matter, inside the orbits 
of Mars or Mercury. We obtain the following expres¬ 
sions for < vj >, on taking into account the orbit 
of Mars (e = 9.34 x IQ-^, i = 1.84, a = 1.52 AU, 
( Horizonsvstem 2015h l: 


< rij >Mars= —7.69s X 10 ^®mas cty ^ 
while for Mercur y (e = 2.05 x 10“^ 


3.87 X 10 ^ AU, ( Horizonsvstemll2015ip : 
< w >Mercury= 3.44s X 10“^®mas cty“^ 


(29) 


= 7.00, a = 


( 30 ) 





















































6 


On using the data obtained by iFienga et al.l (|201ll) . 
AdjMars = —0.04±0.15 mas cty“^; we obt ain s < 1.43x 
10^® kg s“^. As for Mercury, it is (jFienga et al. 


20111) AnjMercury = 0.4 ± 0.6; similarly, we obtain 
s < 2.90 X 1025 kg m2 s-k 

As for planets orbiting inside the ring, we see that 
< rij > in Eq. (E51) is constant, so that it does not 
depend on the orbit of the test particle (even though 
the orbit is not in the plane of the ring): as for its 
magnitude, we obtain 


/ 1AU\^ 

< w >= s 2.9 X 10“^5 / ^ j mas cty 


-1 


(31) 


As a consequence, the spin magntiude s of a hypo¬ 
thetical ring at i? = 1 AU can be co nstrained by 
using the data of Venus, measured by iFienga et al. 
(120nil : Adjvenus = 0.2 ± 1.5 mas cty”^; we obtain 
s < 5.9 X 10^2 kg m2 s“5. Similarly, we can con¬ 
strain the spin of the minor asteroids belt between 
Mars and Jupiter, by considering R = 2.8 AU, and 
using the perihelion of Mars, measured bv IFienga et al. 
(j201lh AdjMars = —0.04 ± 0.15 mas cty“^; we obtain 
s < 8.3 X 10"‘2 kg m2 s“^. 

It is important to explain the meaning of the above 
estimates: indeed, they should be considered just as 
upper limits, useful to evaluate the order of magnitude 
of the effects. In fact, in actual physical situations, the 
GM perturbations due to the ring are present together 
with other effects, such as the Lense-Thirring, the J 2 
effects of the central body and the Newtonian or GE 
effects of the ring; among the latter, there are the tidal 
interactions: in particular, it is possible to show that 
for actual physical situations, the impact of the tidal in¬ 
teractions are much greater than the GM perturbations 
due to ringH Moreover, in the case of planets orbiting 
outside the ring we notice that the expression of the 
GM field (j21l) is the same as the GM field of the central 
body, in terms of its own angular momentum ([8]); in 
particular, the secular variations are t he same as those 
of the classical Lense-Thirring effect (|lorioll200ll) . As 
a consequence, far away from the central body and the 
ring, the total GM field will depend on the sum of the 
angular momenta of the central body and the ring, and 


rough estimate of the ratio of the magnitudes of the GM ac¬ 
celeration 'Vf QM (due to a ring of mass orbiting at distance 
Rr from a central body of mass M) to the tidal acceleration Wy 
(due to a planet of mass rrip orbiting the same central body at dis¬ 
tance Rp) on test particle orbiting at distance r {r < Rr,r < Rp, 

GM mr 

can be written in the form - ^ , which suggests 

/ r 

y Rr Rl 

that the post-Newtonian GM effects are very small. 


it would be very difficult (at least for the chosen ring 
configuration) to set constraints on the angular momen¬ 
tum of the ring. 


5 Conclusions 

In this paper we have focused on the GM field produced 
by rotating rings of matter, orbiting around a central 
body, regarded as a small perturbation of the leading 
gravitational field of the central body. In particular, 
we have considered a thin circular ring, with constant 
matter density, and calculated its field, in the form of 
power law, in the intermediate zone between the cen¬ 
tral body and the ring, and also far away from the ring 
and the central body. Then, we have used the lowest 
order expression of the GM field, both inside and out¬ 
side the ring, to calculate the corresponding perturb¬ 
ing acceleration on the Keplerian orbit of a test parti¬ 
cle, with arbitrary inclination with respect to the ring 
plane, thus extending some previous results. As a pos¬ 
sible application, we have evaluated the impact of the 
GM perturbations on the Keplerian orbital elements, 
to make a comparison with the available data in the 
Solar System: namely, on taking into account the data 
of the planetary ephemeris, we have used the predicted 
perturbations of the orbital elements to give rough esti¬ 
mates on the spin of ring-like structures. These results 
are preliminary: the simple model that we have consid¬ 
ered, in fact, can be used to obtain upper limits on the 
spin of the rings, since we have not taken into account 
the other perturbations that are present. However we 
suggest that, at least in principle, by means of a more 
realistic and systematic analysis of the perturbations, 
it could be possibile to infer more information on the 
spin of ring-like structures. 
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